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ABSTRACT 

Non-urban  were  more  favorable  than  urban  stream  sections  as  habitat 
for  trout  and  held  more  trout.   The  major  habitat  difference  was  amount 
of  instream  solid  overhead  hiding  cover.   Urban  land  modifications  had 
created  unnaturally  straight,  narrow  channels  with  high,  unstable  banks 
with  little  of  the  undercuts  and  woody  debris  that  provide  shelter  for 
fish.   Urban  and  non-urban  sections  did  not  differ  significantly  with 
respect  to  water  velocity,  dissolved  nitrate,  or  amount  of  pools  or 
water  turbulence.   Per  unit  stream  length,  non-urban  sections  averaged 
54%  more  trout  larger  than  20  cm  (8  inches)  and  74%  greater  total  trout 
biomass  than  urban  sections. 

In  both  urban  and  non-urban  areas,  trout  abundance  as  kg/ha  was 
generally  below  the  level  predicted  by  the  Wyoming  Habitat  Quality 
Index  (HQI).   This  could  have  been  due  to  effects  of  angling  or  other 
unmeasured  factors,  to  measurement  errors  or  to  inapplicability  of  the 
HQI  method  to  the  areas  studied.   There  is  evidence  that  altering  the 
HQI  method  to  consider  solid  overhead  hiding  cover  and  pool-turbulence 
hiding  cover  as  separate  variables  rather  than  as  a  total  cover  index 
will  enhance  predictiveness. 

Implications  for  urban  stream  fishery  management  are  discussed. 
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PROJECT  SUMMARY 

Trout  populations  and  habitat  were  analyzed  in  16  urban  and  14 
non-urban  areas  of  four  streams  in  and  near  Bozeman  and  Livingston, 
Montana.   The  non-urban  stream  sections  had  generally  undergone  less 
artificial  alteration. 

Non-urban  stream  on  average  held  74%  greater  weight  (biomass)  of 
trout  per  unit  length  of  channel  than  did  urban  stream.   Number  of 
over-20-cm  (over-8-inch)  trout  per  unit  channel  length  was  54%  greater 
in  non-urban  than  in  urban  areas. 

The  urban  parts  of  streams  were  clearly  less  favorable  habitat  for 
trout  in  certain  key  respects.   The  most  striking  habitat  difference 
detected  between  urban  and  non-urban  areas  was  that  urban  stream  had 
significantly  less  hiding  cover  for  trout — hiding  cover  measured  as 
amount  of  solid  material  in  the  water  or  within  a  few  centimeters  above 
it  that  could  provide  overhead  concealment  for  fish.   Also,  in  urban 
stream,  channel  width  was  less  and  amount  of  eroding  bank  greater.   It 
was  apparent  that  urban  landfill  along  stream  banks  had  created  straight, 
narrow  channel  with  high,  unstable  banks.   Such  channels  tended  to  lack 
undercut  banks  that  could  provide  shelter  for  trout.   Within  such 
straightened  and  constricted  waterways,  the  greater  forces  during  high 
water  may  have  swept  away  logs  and  other  woody  debris  that  would  have 
formed  trout  cover.   Also,  people  may  have  removed  debris  to  tidy  the 
appearance  of  urban  channels  or  to  help  them  conduct  flood  water  more 
rapidly,  thereby  destroying  trout  habitat. 
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Urban  and  non-urban  stream  did  not  differ  significantly  with  respect 
to  water  velocity,  nitrate  content,  or  amount  of  cover  that  was  in  the  form 
of  pools  or  water-surface  trubulence. 

Trout  abundance  expressed  as  kilograms  per  hectare  was  generally  below 
the  level  indicated  as  potential  according  to  the  Wyoming  Habitat  Quality 
Index  (HQI)  in  both  urban  and  non-urban  areas.   This,  as  well  as  poor 
coorelation  betwen  HQI  values  and  measured  trout  abundance  may  have  been 
due  to  various  possible  circumstances:   (1)  suppression  of  trout  abundance 
by  some  unmeasured  factor  or  factors,  such  as  angling,  pollution  or 
artificial  diminution  of  streamflow  discharge,  (2)  inapplicability  of  the 
HQI  method  to  the  kinds  of  stream  in  the  study,  (3)  error  in  estimation 
of  habitat  and/or  trout  population  variables,  or  (4)  inappropriateness  of 
biomass  per  unit  stream  surface  area  as  a  measure  of  trout  abundance. 

A  second  type  of  regression  analysis  of  association  between  trout 
abundance  and  habitat  variables  indicated  that  the  lumping  of  all  cover 
types  (pool,  turbulence  and  overhead)  which  is  involved  in  the  HQI  method 
may  have  been  a  major  source  of  poor  predictiveness.   Multiple  correlation 
was  far  higher  when  overhead  cover  and  pool-and-turbulence  cover  were 
treated  as  separate  variables.   This  analysis  indicated  that  71%  of 
variation  in  number  of  over-20-cm  trout/km  was  attributable  to  the 
combined  variation  in  the  following  habitat  factors  in  this  order  of 
importance:   (1)  mean  water  velocity  in  stream  reach,  (2)  ratio  of  late 
summer  streamflow  discharge  to  mean  annual  discharge,  (3)  amount  of 
nolid  overhead  cover,  ('I)  maximum  summer  water  temperature,  and  C?) 
amount  of  aquatic  vegetation  in  the  stream. 


3 

Trout  abundance  in  some  of  the  relatively  unaltered  parts  of  Bozeman 
Creek  compared  favorably  with  that  in  the  more  densely  populated  small 
trout  streams  of  the  same  geographic  area. 

The  results  suggest  that  trout  abundance  can  be  maintained  in  streams 
flowing  through  areas  under  urban  development  if  the  natural  form  and 
vegetation  of  the  stream  banks  and  stream  bed  are  not  altered  in  the  ways 
commonly  associated  with  urbanization.   If  urban  changes  in  land  form 
are  kept  well  away  from  the  immediate  riparian  area,  the  natural  channel 
shape  and  its  natural  accumulations  of  living  and  dead  vegetation  will 
furnish  cover  for  substantial  populations  of  trout.   It  is  especially 
important  not  to  straighten  channels ,  not  to  remove  certain  kinds  of  bank 
vegetation  (such  as  high  grass  and  low  brush),  and  not  to  conduct  excessive 
removal  of  downed  logs  and  other  woody  debris  from  stream  channels. 

Much  could  be  done  to  increase  trout  abundance  in  physically  damaged 
urban  streams  by  restoring  channel  form  and  vegetation  to  resemble  the 
natural  situation.   Even  in  parts  of  streams  that  remain  unaltered,  habitat 
may  often  be  enhanced  and  trout  populations  increased  by  creating  more 
instream  cover  for  trout  than  presently  exists. 
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INTRODUCTION 

Expanding  urbanization  and  its  effects  of  unfavorably  reshaping  the 
channels  of  urban  streams  is  a  serious  fisheries  problem  that  has  long 
seemed  obvious  to  biologists  but  may  not  be  as  apparent  to  others.   Urban 
stream  alterations  typically  regarded  as  detrimental  to  fish  habitat 
include  channel  straightening  and  other  unnatural  relocation;  excessive 
stream  widening  or  channel  constriction  through  landfill,  bridges  and 
culverts;  impoundments;  elimination  of  pools,  riffles  and  biologically 
productive  side  channels;  removal  of  instream  cover  used  by  fish;  and 
construction  of  unnatural  structures  such  as  bulkheads,  walls  and  certain 
kinds  of  deflectors  in  and  along  streams.  This  situation  has  not  been 
thoroughly  studied  in  quantitative  terms,  in  large  part  perhaps  because 
it  has  seemed  so  blatantly  obvious  to  specialists.   Yet,  when  it  comes 
to  making  policy  decisions  on  urban  development  along  streams,  quantitative 
data  are  needed. 

Determining  components  of  fish  habitat  that  are  adversely  affected  by 
practices  similar  to  these,  although  not  necessarily  classifiable  as  urban, 
has  been  addressed  in  several  studied.   Whitney  and  Bailey  (1959)  and  Elser 
(1967)  both  found  that  channel  alterations  involved  in  highway  construction 
caused  significant  decrease  in  the  salmonid  carrying  capacity  of  two  Montana 
streams.   In  particular,  Elser  determined  the  amount  of  cover  per  unit  area 
of  stream  to  be  substantially  less  and  the  occurrence  of  areas  with  more 
shallow,  fast  reaches  greater  in  altered  stream  sections.   In  a  stream- 
straightening  project  in  Rocky  Creek,  Gallatin  County,  Montana,  changes 
in  channel  morphometry  and  loss  of  in-streara  and  bank  cover  were  responsible 


5 

for  reduced  trout  abundance  (Wells  1977).   Still  other  studies  have  found 
trout  populations  to  be  limited  by  quality  of  physical  habitat  (Boussu 
1954,  Kalleberg  1958,  Lewis  1969,  and  Newman  1956). 

Our  study  was  intended  to  evaluate  habitat  quality  by  comparing 
measurements  of  various  habitat  attributes  between  urbanized  and  non-urbanized 
sections  of  stream  and  analysing  them  for  correlation  with  abundance  of 
trout.   While  it  is  almost  impossible  to  locate  pristine,  unaltered  streams 
in  and  near  urban  areas,  one  can  probably  make  adequate  analyses  by  comparing 
altered  and  "less  altered"  parts  of  streams  in  and  near  towns  within  generally 
urbanized  areas. 

In  urban  creeks  of  the  Bozeman  area  in  Montana,  stream  problems  have  been 
studied  from  water  quality/pollutional  standpoints  (Anderson  1977;  Blue 
Ribbons  of  the  Big  Sky  Country  APO  1979).   The  emphases  of  these  studies 
were  on  sediment  pollution  and  on  pollutants  hazardous  to  human  health. 
The  urban  reaches  of  streams  were  not  analysed  with  respect  to  physical 
suitability  as  habitat  for  fish. 

The  hypotheses  selected  for  this  study  were  (1)  that  key  habitat 
characteristics  are  less  favorable  to  trout  in  altered  than  in  less-altered 
parts  of  urban-area  streams  and  that  (2)  trout  are  less  abundant   in  the 
more  altered  parts  of  urban-area  streams. 

DESCRIPTION  OF  STUDY  AREAS 
Thirty  channel  reaches  (stations)  totalling  3,772  meters  of  stream  were 
analysed  for  trout  habitat  attributes  and  trout  abundance  in  four  creeks  that 
have  undergone  varying  degrees  of  artificial  alteration  in  the  course  of  urban- 
ization (Table  1).   Three  of  the  streams  flow  through  the  city  of  Bozeman  (Fig.  1) 


Table  1.   Number  of  study  stations,  total  length  of  stream  studied,  and 
trout  species  present  in  the  study  streams. 
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30 
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Figure  1.   Streams  in  the  vicinity  of  Bozeman,  Montana. 


Figure  2.   F'' eshtnan  Creek  and  Livingston,  Montana. 
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in  Gallatin  County:   Bozeman  Creek  (18  stations  of  ca.  100-200  m  each,  distributed 
over  6.4  km  of  stream),  Mathew  Bird  Creek  (7  stations)  and  Figgins  Creek  (2 
stations).   One  stream,  Fleshman  Creek  (3  stations)  flows  through  the  city  of 
Livingston  (Fig.  2)  in  Park  County.   Details  of  station  locations  are  in  Appendix  I 

Urban  stations  did  not  include  the  most  drastically  altered  parts  of  some 
streams — those  parts  that  had  been  completely  enclosed  in  culvert  for  one  or 
more  city  blocks.   It  would  have  been  too  dangerous  to  electrofish  and  perhaps 
also  to  do  other  sampling  in  these  sections. 

Study  station  lengths  ranged  from  89  to  201  m,  except  for  one  that  was 
300  m  long.   Seventy-two  per  cent  of  stations  were  between  90  and  100  m  long. 

Bozeman  (Sourdough)  Creek  originates  at  the  outlet  of  Mystic  Reservoir 
(elevation  1950  m)  in  the  north  end  of  the  Gallatin  Range  in  southwestern 
Montana  and  flows  northwest  for  26.6  km,  dropping  510  m  (19.2  m/km)  to  its 

mouth  (elevation  1440  m)  at  the  East  Gallatin  River  below  Bozeman.   The  Bozeman 

2 
Creek  drainage  basin  covers  about  168  km  .   Average  annual  precipitation  for 

the  basin  is  6l  cm  (24  inches). 

Mathew  Bird  Creek  (Spring  Creek  on  the  Bozeman  topographic  quadrangle  of 
1953)  is  tributary  to  Bozeman  Creek  at  about  stream  km  3-5  at  elevation 
1475  m.   Elevation  of  the  source  is  1585  m.  Stream  length  is  about  3.3  km. 

Figgins  Creek  (Middle  Creek  Ditch),  enters  Mathew  Bird  Creek  at  stream 
kilometer  1.6.   It  originates  from  Hyalite  Creek  in  T3S,R8E,Sec.  3- 

Fleshman  Creek  heads  in  eastern  foothills  of  the  Bridger  Range  in 
Park  County  and  flows  southeast  for  14.8  km  before  discharging  into  the 
Yellowstone  River  in  the  city  of  Livingston.   Its  elevations  at  origin  and 
mouth  are  1707  m  and  1372  m,  respectively.   The  average  annual  precipitation 
for  the  Fleshman  Creek  drainage  basin  is  40.6  cm  (16  inches). 
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Large  parts  of  these  four  streams  are  bordered  by  agricultural,  industrial 
and  municipal  lands.   Land  uses  include  ranching,  small  grain  production,  small 
industry,  and  urban  development — and  in  the  headwaters  of  Bozeman  and  Fleshman 
Creeks,  forest  uses.   A  shift  from  agricultural  to  urban-municipal  use  is 
rapidly  occurring  in  many  parts  of  the  bordering  lands,  with  increase  in 
artificial  alterations  of  stream  channels  and  banks. 

The  range  of  approximate  annual  mean  discharges  over  all  study  sites  was  0.04 
m  /sec  to  1.46  m  /sec.   Average  critical  period  discharge  (late  summer  flow  from 
August  1  to  September  15)  ranged  from  0.034  m  /sec  to  0.80  m  /sec. 

METHODS 
Selection  of  Study  Stations 

Potential  study  streams  were  suggested  by  fishery  biologists  of  the 
Montana  Department  of  Fish,  Wildlife  and  Parks.   Identification  of  deliberate 
stream  course  relocations  done  in  recent  decades  was  accomplished  by  inspection 
of  aerial  photos  (scale:   12  inches/mile)  taken  in  1937,  1954  and  1977. 

Each  stream  was  divided  into  reference-station  reaches  by  measuring 
and  marking  of  100-m  reaches  contiguously  from  the  mouth  upstream  for  several 
kilometers  including  both  urban  and  non-urban  areas,  but  with  some  gaps  in  the 
system  of  stations,  owing  to  access  problems.   Within  each  stream,  the  stations 
were  classified  as  being  either  urban  (definitely  within  an  urban-impacted  area) 
or  control  (non-urban — although  perhaps  "suburban"  in  some  cases — and  relatively 
less  impacted  by  artificial  alteration).   Within  each  category,  a  sample  of 
stations  was  selected  for  study. 

On  Bozeman  Creek,  many  of  these  stations  were  contiguous  in  groups  of  two 
or  three.   For  some  analyses,  contiguous  stations  were  treated  as  combined 
areas,  called  "sections,"  with  the  data  lumped  within  sections. 
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Estimates  of  Trout  Population  and  Biomass 

Between  August  1  and  September  8,  1982,  electrofishing  for  mark-recapture 
estimates  of  trout  populations  was  done  in  each  of  the  30  stations  selected 
for  study.   The  electrofishing  unit  consisted  of  a  240-volt,  3500-watt  alternator 
and  Coffelt  control  box  (150  to  600  volts  DC)  mounted  in  a  canoe  and  towed 
upstream  by  a  crew  using  two  positive  electrodes  ahead  of  the  unit,  a  negative 
electrode  trailing  in  the  water  behind  the  canoe — except  in  very  shallow 
stations,  where  long-line  gear  was  used  instead  of  the  canoe-mounted  unit. 
Collected  fish  were  anesthetized  in  MS-222  (Tricane  Methanesulfonate ) ,  weighed 
to  the  nearest  2  grams,  and  measured  to  the  nearest  1  mm  in  length.   The  lower 
caudal  fin  was  clipped;  when  two  300-meter  sections  were  contingous,  the  upper 
caudal  fin  of  the  upstream  section  trout  was  clipped.   Fish  were  carried  in 
large  frame  nets  downstream  to  the  bottom  of  the  station  from  which  they  had 
been  caught  and  released.  Two  weeks  were  allowed  between  the  marking  and 
recapture  runs. 

The  population  estimate  for  each  station  (and  3-station  section)  was 
calculated  by  the  modified  Petersen  mark-recapture  formula  of  Ricker  (1975). 
Fish  were  grouped  into  20-mm  size  classes,  100  to  119   mm  being  the  smallest. 
Where  numbers  were  too  low  to  calculate  a  reliable  pouplation  estimate  for  each 
species  within  a  station,  data  were  lumped  either  by  species  or  by  combining 
several  stations,  and  the  estimates  calculated  then  reapportioned  according 
to  total  of  fish  marked  on  first  run  plus  the  unmarked  (new)  fish  caught  on  the 
second  run. 

Biomass  was  calculated  by  multiplying  length-group  estimates  of  trout 
numbers  by  mean  weight  for  the  group.   Mean  weight  of  length  group  was 
determined  from  length-weight  regression  equations  for  each  species  in  each 
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station.   Standing  crop  in  kilograms  and  numerical  density  of  over-20-cm 
trout  were  calculated  on  per-hectare  and  per-stream-kilometer  bases. 

The  95%  confidence  intervals  (Ricker,  1975:  8l)  for  the  estimates  of 
numbers  of  trout  over  20  cm  were  also  calculated  for  each  station  and  section 
and  converted  to  density. 
Habitat  Measurements 

Stream  habitat  attributes  were  measured  during  the  low  flow  period  of 
September  1982  to  March  I983,  following  the  prcedures  described  by  Binns 
(1979,  1982)  for  the  Wyoming  Habitat  Quality  Index.   Habitat  variables  measured 
or  evaluated  were  maximum  summer  water  temperature,  stream  discharge,  channel 
width,  percent  of  eroding  banks,  length  of  thalweg,  mid-channel  length,  amount 
of  submerged  aquatic  vegetation  (substrate),  nitrate  nitrogen  level,  mean 
water  velocity  through  station,  and  percent  cover.   The  cover  designation 
included  any  of  the  following  features  that  occurred  in  water  at  least  15  cm 
deep:  undercut  banks,  instream  and  closely  overhanging  terrestrial  vegetation, 
instream  debris  (brush,  logs  and  "snags"),  pockets  of  surface  turbulence,  and 
pools.   Pools  were  identified  as  abrupt  increases  in  water  depth.   The  length 
and  width  of  these  cover  features  were  measured,  summed  and  converted  to  per 
cent  of  total  stream  surface  area  in  the  station.   Subsequently,  the  categories 
of  pool  and  surface  turbulence  pockets  were,  for  more  detailed  analyses, 
separated  from  the  other  cover  types — which  could  be  categorized  as  "solid 
overhead  cover." 

Mean  water  velocity  through  station  was  measured  by  timing  an  injection  of 
rhodamine  dye  (Binns  1982). 

Nitrate  nitrogen  levels  were  measured  by  the  ultraviolet  spectrophotometric 
method  as  described  in  Standard  Methods  for  Examination  of  Water  and  Wastewater 
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(1976). 

Annual  stream  flow  variation  was  determined  by  the  use  of  past  discharge 
records  and  estimates  from  governmental  agencies .   Late  summer  streamf low 
discharge  for  each  stream  was  measured  with  a  gurley  "pigmy"  flow  meter. 
Stream  discharge  gauging  records  were  available  only  for  Bozeman  Creek.   The 
USDA  Soil  Conservation  Service  operated  a  gauge  at  stream  kilometer  5.6  in  1977 
and  1978. 
Data  Analysis 

The  mean  values  of  urban-  and  control-section  data  on  habitat  and  trout 
population  variables  were  compared.   The  data  were  examined  for  significant 
differences  via  the  non-parametric  Mann-Whitney  U  test. 

Data  on  habitat  variables  were  entered  into  the  Wyoming  Habitat  Quality 
Index  Model  II  (Binns  1979).   The  resultant  predictions  of  trout  abundance  in 
kilograms  per  hectare  were  compared  with  the  actual  measurements  of  trout 
abundance. 

The  relationship  of  trout  abundance  to  habitat  variables  was  further 
analysed  by  multiple-regression  modelling.   The  logarithms  (base  10)  of  habitat 
data  were  entered  as  independent  variables.  The  logarithms  of  four  expressions 
of  trout  abundance  (number  of  over-20-cm  trout  per  stream  kilometer  and  per 
hectare  and  kilograms  of  all  trout  per  stream  kilometer  and  per  hectare)  were 
entered  separately  as  dependent  variables.   In  the  step-wise  procedure,  the 
F- value  set  for  exclusion  of  independent  variables  from  the  model  was  1.0. 

RESULTS  AND  DISCUSSION 

Comparisons  between  Urban  and  Control  Areas 

The  most  striking  difference  detected  between  urban  and  control  areas  of 
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the  streams  was  in  terms  of  amount  of  solid  overhead  concealment  cover  for 
trout — cover  submerged  in  the  stream  or  within  a  few  centimeters  above  the 
water  surface.   Averaged  over  all  30  study  stations,  control  areas  had  60%  more 
such  cover  than  did  urban  areas — 8.33%  vs  5.22%,  a  difference  significant  at 
the  97.5%  confidence  level  (Table  2).   For  Bozeman  Creek  alone,  the  difference 
was  even  more  pronounced.   Urban  areas  contained  especially  low  amounts  of  cover, 
slightly  less  than  1%,  and  the  control  areas  did  not  have  particularly  high 
amounts,  5.3%,  but  this  was  a  several-fold  difference  or  485%  more  in  control 
areas,  a  difference  significant  at  the  99.8%  level. 

It  can  be  inferred  that  reduction  of  instream  concealment  cover  is  a  major 
impact  of  urbanization  along  trout  streams  of  the  types  included  in  our  study. 
Wells  (1977)  also  found  that  cover  reduction  was  strongly  associated  with 
channel  straightening  in  another  stream  near  Bozeman,  although  that  was  not  a 
consequence  of  urbanization. 

The  only  other  measured  habitat  variables  having  significant  urban-vs-control 
differences  at  the  90%  confidence  level  or  higher  were  channel  width,  which 
was  narrower  in  urban  areas  (96.3%  confidence  level),  and  amount  of  eroding 
bank — 34%  in  urban  areas  vs  23%  in  control  areas,  a  difference  of  48%  relative 
to  the  control  figure  and  significant  at  the  90%  confidence  level  (Table  2). 

Part  of  the  difference  in  channel  width  could  be  an  artifact  of  the 
longitudinal  distribution  of  study  stations  within  the  stream  systems.   This 
possible  effect  has  not  yet  been  analysed.   However,  we  suspect  that  encroach- 
ment on  the  channel  by  land  fill  involved  in  urban  development  may  have  been  a 
major  influence  in  narrowing  the  urban  stream  sections. 

The  higher  degree  of  bank  erosion  in  urban  areas  may  represent  decreased 
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Table  2.   Unweighted  mean  values  for  habitat  and  trout  population  variables  in 
urban  (U)  and  control  (C)  stations  of  study  streams.   Values  of  p  are 
from  Mann-Whitney  U  tests. 


Variable 


Station 
type 


Bozeman 

Creek 
(n=9U,9C) 


Mathew 

Bird 
Creek 
(n=4U,3C) 


Figgins     Fleshman 
Creek       Creek 
(n=lU,lC)    (n=2U,lC) 


Habitat   Variables 


Channel       U 
width  (m)      C 


4.84 

7.01 

(p<.001) 


2.31 
2.34 
(ns) 


1.40 
1.58 


2.15 
1.73 


3.65 

5.24 

(p=.047) 


Solid  over-    U 
head  cover  (%) 


0.91 

5.32 

(p<.002) 


9.24 

17.53 

(ns) 


30.42 
13.33 


3.94 
2.84 


5.22 

8.33 

{p=.025) 


Pool,  turbu-   U 
lence  cover (%)  C 


7.13 
4.78 
(ns) 


8.52 
7-21 
(ns) 


0.50 
5.47 


7.38 
6.38 


7.09 

5.46 

(p=.377) 


Eroding 
banks    (%) 


32.85 

19.07 

(ns) 


24.01 
1.95 
(ns) 


12.65 
9.00 


70.3 
126.0 


33.8 
22.6 
(p=.101) 


Water 

velocity 

(cm/sec) 


79.5 
77.2 
(ns) 


43.5 
22.9 
(p<.02) 


12.2 
12.2 


32.6 
26.0 


60.4 
57.3 
(p=.790) 


NO   -N 
(m|/l) 


0.420 
0.371 
(ns) 


0.612 
0.338 


0.325 
0.338 


0.414 
0.429 


0.393 

0.431 

(p=.822) 


Trout      Variables 


Number  of   Trout   20   cm  and   Larger 


Per  stream  km       U  323  188  64 

C  523  79  241 

(p<.02)  (ns) 


244 
376 
(p=.154) 


Per  stream  ha   U     673        827         456 
C     754         346        1524 
(ns)        (ns) 


421 
483 


666 
702 
(p=.984) 


Kilograms  of  All  Trout 


Per  stream  km 

U 

56.9 

35.4 

15.3 

20.1 

44.3 

C 

104.1 
(p<.01) 

24.7 
(ns) 

51.0 

18.6 

77.2 
(p=.052) 

Per  stream  ha 

U 

118.5 

155.1 

109.5 

93.1 

123.9 

C 

152.4 
(ns) 

107.5 
(ns) 

322.9 

107.5 

151.8 
(p=.355) 
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channel  stability,  a  hydraulic  response  of  the  streams  to  construction  of 

unnaturally  straight  and  narrow  channels  with  high,  steep  and  constricting 

banks.   There  had  apparently  also  been  some  past  artificial  channel  straightenings 

in  some  of  the  control  stations.   Had  that  not  been  the  case,  perhaps  the  difference 

in  amounts  of  bank  erosion  between  urban  and  control  areas  would  have  been  even 

greater. 

Most  of  the  erosion  was  at  high  elevation  on  banks,  apparently  coinciding  with 
spring  flood  stage.   Urban  station  4  on  Mathew  Bird  Creek,  while  having  only 
9.76%  erosion  of  banks,  showed  obvious  potential  for  erosion  due  to  urban 
development.   With  construction  of  condominiums  occurring  within  30  meters 
or  less  of  the  stream,  all  of  the  riparian  willows  and  high  grasses  had  been 
removed.   Where  natural  deeply  undercut  bank  cover  had  been  formed  by  the 
stream  action,  major  portions  of  it  were  beginning  to  slump  into  the  stream, 
leaving  bare  soil  exposed.   Fleshman  Creek's  stations  (including  the  control 
station  which  had  apparently  been  ditched  some  years  ago)  had  the  highest 
average  percent  of  eroding  banks  (88.89%).   We  do  not  believe  eroding  stream 
margins  to  be  of  direct  detriment  to  trout  in  the  immediate  vicinity  of  the 
erosion  but  suspect  they  may  be  indicative  of  unfavorable  action  of  currents  at 
high  flow  in  unnatural  channel  configurations ,  and  there  may  often  be  damage 
to  downstream  areas  via  siltation  of  spawning  and  food-producing  stream  bed. 

The  urban  areas  had  somewhat  greater  measured  amounts  of  stream  area  in 
pools  or  having  the  water  surface  broken  by  turbulence  that  might  provide 
concealment  cover  for  fish.   However  the  urban-vs-control  differences  in  this 
variable  were  not  statistically  significant  at  the  95%  level  of  confidence. 
There  was  even  less  difference  in  the  water  velocity  and  nitrate-nitrogen 
variables. 
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Trout  abundance  was  higher  overall  in  control  stations  than  in  urban  stations 
with  respect  to  all  four  expressions  of  abundance  that  were  calculated  (Table  2).   | 
Trout  abundances  expressed  per  unit  of  stream  length  showed  much  greater  statistical  I 
significance  of  urban-vs-control  differences  than  did  per-unit-area  expressions. 
In  terms  of  kilograms  per  stream  kilometer,  control  areas  had  lk%   more  trout  than 
urban  areas  did,  a  difference  significant  at  the  97-5%  confidence  level.   In 
Bozeman  Creek,  the  control  sections  averaged  83%  more  kg  of  trout  per  km,  a 
difference  significant  at  the  99%  level.   The  urban-vs-control  differences 
in  terms  of  number  of  over-20-cm  (over-8-inch )  trout  per  kilometer  were  significant 
at  the  85%  level  overall  and  at  the  98%  level  for  Bozeman  Creek,  which  had  the 
majority  of  stations. 

It  can  be  concluded  that  the  results  are  consistent  with  the  hypotheses 
of  this  study:   (1)  that  key  habitat  variables  are  less  favorable  for  trout  in 
altered  (urban)  than  in  less  altered  stream  sections,  and  (2)  that  trout  are 
less  abundant  in  the  more  altered  sections.   The  following  sections  of  this 
report  further  support  these  conclusions  and  provide  some  information  toward 
identifying  possible  causative  processes. 
Application  of  the  Wyoming  HQI 

Application  of  the  Wyoming  Habitat  Quality  Index  Model  II  (Binns  1979)  to 
the  data  from  all  30  study  stations  yielded  poor  correlation  (r  =  0.228)  between 
predicted  and  actual  standing  crops  of  trout  expressed  as  kg/ha  (Table  -3,  Figure 
3),  however  it  appears  that  if  data  from  study  sections  (combinations  of  contiguous 
stations)  were  plotted  instead  of  stations,  variability  would  be  reduced  and  the 
fit  would  be  somewhat  tighter.   (This  analysis  remains  to  be  done.)   It  would  still 
be  the  case,  however,  that  in  the  great  majority  of  cases,  predicted  values  fall 
below  the  actual  values  of  standing  crop. 
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Table  3.   Comparison  of  trout  standing  crop  predictions  by  the  Wyoming  HQI 
method  (Model  II)  against  measured  values. 

~~~  Standing  Crop  of  Trout  (kg/ha )~ 

Station  Type*  Predicted  Measured 

Bozmean  Creek 

141 
91.4 
9  C  253  175 

13  C  229  306 

14  C  237  117 

15  C  229  119 

19  U  26.2  83.2 

20  U  253  181 

21  U  179  81 

22  U  148  ■  123 

23  U  148  216 

24  U  226  125 

26  U  116  97.1 

27  U  148  93.9 

28  U  148  73-0 
36  C  148  227 

60  C  214  95.7 

61  C  296  100 


1  U  307  160 

2  U  181  179 

3  U  285  148 

4  U  148  133 

5  C  265  34. 

6  C  265  161 

7  C  297  127 


1  U  259  109 

2  C  223  323 


1  U  61.2  96.2 

2  U  61.2  90.0 

3  C  61.2  108 


c 

284 

c 

150 

c 

253 

c 

229 

c 

237 

c 

229 

u 

26.2 

u 

253 

u 

179 

u 

148 

u 

148 

u 

226 

u 

116 

u 

148 

u 

148 

c 

148 

c 

214 

c 

296 

Mathew  Bird  Creek 

u 

307 

u 

181 

u 

285 

u 

148 

c 

265 

c 

265 

c 

297 

Figgins  Creek 

u 

259 

c 

223 

Fleshman  Creek 
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Figure  3.  Relationship  of  standing  crops  of  trout  as  measured  in  the 

study  stations  and  as  predicted  by  the  Wyoming  Habitat  Quality 
Index  Model  II. 
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This  could  be  interpreted  in  several  ways.   The  HQI  model  may  not  be  applicable 
to  the  trout/habitat  relationships  in  the  kinds  of  streams  included  in  our  study. 
The  poor  fit  of  the  data  might  also  be  due  to  deficient  accuracy  or  precision 
in  our  measurement  of  some  habitat  variables  or  of  standing  crop  of  trout.   Some 
of  the  habitat  measurements  called  for  in  the  HQI  method  seem  more  subjective  than 
would  be  desirable  (e.g.,  vegetative  substrate)  or  anotherwise  difficult  to 
accurately  determine  (e.g.  annual  streamflow  variation),  but  this  may  be 
compensated  for  at  least  partially  by  the  transformation  of  measurements  to  class 
ratings  involved  in  the  procedure  (Binus  1979). 

Comparison  of  the  results  of  the  HQI  method  against  results  of  the  multiple- 
regression  analysis  in  the  following  section  of  the  report  would  seem  to  indicate 
that  much  of  the  unpredictiveness  of  the  HQI  in  our  streams  lay  in  the  procedure 
of  combining  all  forms  of  cover,  i.e.,  failure  to  distinguish  between  pool-and- 
turbulence  cover  and  solid  overhead  cover.   Comparison  of  the  predictiveness  of 
the  HQI  method  in  our  streams  against  the  better  predictiveness  of  the  multiple- 
regression  analyses  desribed  below  may  also  point  to  an  inappropriateness  of  using 
biomass  per  unit  of  stream  surface  area  as  an  expression  of  trout  abundance. 

It  is  far  more  important  to  consider  the  possible  inference  from  the 
lower-than-predicted  distribution  of  actual  trout  abundances  that  some 
unmeasured  influence  is  preventing  trout  in  these  streams  form  being  as  abundant 
as  the  HQI  predicts  they  should  be,  i.e.,  as  the  physical  habitat  would  allow. 
Actual  abundance  is  lower  than  HQI-predicted  abundance  particularly  in  the 
stations  with  highest  predicted  values,  i.e.,  highest  habitat  rating  (Fig.  3). 
Candidate  variables  for  such  abundance-depressing  effect  would  be  water  pollution 
and  intense  angling  harvest.   Neither  were  measured  in  this  study,  however,  urban 
water  pollution  problems  were  revealed  in  a  previous  study  on  Bozeman  Creek  (Blue 
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Ribbons  of  the  Big  Sky  Country  APO  1979). 

It  could  also  be  that  streamflow  variability  was  not  adequately  estimated. 
The  summer  low  flow  may  often  be  less  due  to  irrigation  withdrawal  than  it  was 
when  measured  in  1982,  a  fairly  wet  summer.   Also,  there  may  be  withdrawals 
for  municipal  water  supply,  making  winter  low  flow  severe  in  the  Bozeman  Creek 
system.   This  was  not  analysed  in  our  study. 

Insufficient  reproduction  of  trout  could  be  another  factor  preventing 
standing  crop  from  reaching  the  potential  indicated  by  the  HQI.   Although 
age-structure  analyses  of  the  trout  populations  have  not  yet  been  done,  the  size 
structures  of  the  Bozeman  Creek  rainbow  trout  population,  the  main  fish  in  that 
stream,  appear  to  indicate  fewer  fish  of  age  I  and  II  than  of  age  III  or  IV.   This 
is  the  reverse  of  the  situation  that  must  exist  in  a  population  being  replenished 
by  local  reproduction.   It  is  likely  that  the  rainbow  trout  population  of  our 
study  area  on  Bozeman  Creek  consists  largely  of  immigrants  from  upstream  or  from 
the  East  Gallatin  River,  downstream.   The  combination  of  immigration,  body  growth 
and  whatever  reproduction  exists  within  the  study  may  not  be  creating  enough 
biomass  to  saturate  the  habitat.   If  low  reproductive  rate  is  also  an  influence 
for  brook  trout  in  Bozeman  Creek  and  in  the  other  study  streams  (it  is  the 
predominant  fish  in  the  other  streams),  it  is  not  as  strongly  the  case.   Size 
distributions  of  brook  trout  in  most  study  stations  indicate  that  a  more  normal 
age  structure  probably  exists. 
Multiple-regression  Modelling  of  Trout  and  Habitat  Relationships 

When  logarithms  of  the  habitat  variables  involved  in  the  HQI  (but  not        ~^ 
transformed  to  class  ratings)  were  entered  as  independent  variables  into  stepwise 
multiple  regression  against  logarithms  of  each  of  the  four  trout  abundance  expression; 
used  separately  as  dependent  variables,  stronger  correlations  resulted  than  in  the 
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HQI  model.   Further,  when  the  HQI  cover  variable  was  separated  into  a  solid-overhead- 
cover  component  and  a  pool-and-turbulence-cover  component  and  these  new  components 
entered  as  independent  variables  along  with  the  rest,  much  stronger  correlation 
yet  was  obtained . 

In  a  set  of  stepwise  multiple  regressions  involving  (1)  all  30  study  stations, 
(2)  control  stations  only,  and  (3)  urban  stations  only,  mean  channel  width  was 
the  first  variable  entered,  based  on  its  high  initial  influence  relative  to  that 
of  other  independent  variables.   However,  as  other  variables  were  entered,  mean 
width  rapidly  lost  significance,  its  F-value  falling  below  the  predetermined 
rejection  level  of  1.0.   It  is  inferred  that  channel  width  was  meaningless  in 
describing  the  effect  of  habitat  on  trout  abundance.   Another  stepwise  multiple 
regression  was  run  with  the  channel-width  variable  omitted,  and  the  correlation 
improved  (Table  4). 

Number  of  over-20-cm  trout/km  was  consistently  the  dependent  variable  having 
tighest  correlation  (Table  4).   We  infer  that  the  habitat  variables  are  related 
to  trout  abundance  in  a  linear,  rather  than  areal  fashion.   Trout  abundance  per 
unit  stream  surface  area  is  likely  to  be  poorly  correlated  with  habitat  quality 
because  trout  orient  strongly  to  instream  cover  (Hunt  1971,  Wesche  1976, 
Devore  and  White  1978,  Enk  1977)  which  tends  to  be  concentrated  along  channel 
margins,  hence  is  a  rather  linear  variable. 

With  the  channel-width  variable  omitted  from  multiple  regression,  the 
dependent  trout  abundance  variable  most  highly  correlated  with  habitat  variables 

was  still  the  number  of  over-20-cm  trout/km  (Table  4).   For  the  model  involving 

2 
all  30  stations,  the  adjusted  r  values*  indicated  that  variation  in  habitat 


2    2  2 

♦Adjusted  r  =  r  -  [ (K-1) / (N-k) ][l-r  ],  where  K  =  number  of  independent  variables 

and  N  -   number  of  cases. 
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24 
variables  accounted  for  71%  of  variation  in  number  of  trout/km.   The  habitat 
variables  in  that  model  were  water  velocity,  late  summer  flow/mean  flow, 
solid  overhead  cover,  max  summer  water  temperature,  and  vegetative  substrate 
ratings  (Appendix  III). 

For  the  multiple  regression  models  involving  urban  or  control  stations  only, 
the  r  values  were  even  higher  and  the  solid-overhead-cover  variable  dropped  out 
(Appendix  III),  as  it  had  been  the  most  significantly  differing  variable  between 
the  urban  and  control  categories  of  stream  (Table  2).   Percent  pool-and-turbulence 
cover  and  percent  eroding  banks  entered  only  into  the  urban  model.   Water  velocity 
was  the  first  variable  to  enter  into  all  three  models  (Appendix  III)  after  channel 
width  was  omitted. 
Other  Observations  on  Fish  Populations 

The  species  of  fish  collected  in  the  study  sites  were  rainbow  trout  (Salmo 
gairdneri ) ,  brook  trout  (Salvelinus  fontinalis ) ,  brown  trout  (Salmo  trutta) , 
longnose  sucker  (Catostomus  catostomus) ,  white  sucker  (Catostomus  commersoni) , 
mountain  sucker  (Catostomus  platyrhynchus)  ,  and  mottled  sculpin  (Cottus  bairdi ) . 

Standing  crop  per  unit  stream  area  ranged  from  34  kg/ha  in  station  5  of 
Mathew  Bird  Creek  to  a  high  of  305  kg/ha  in  station  13  of  Bozeman  Creek.  In 
terms  of  numbers  of  over-20-cm  trout  per  kilometer,  abundance  ranged  from  30 
to  611  in  the  same  stations  (Table  5).  Grouping  of  the  stations  into  longer 
stream  sections  reduced  the  confidence  intervals  of  trout  population  estimates 
(Table  6),  enabling  clearer  indication  of  relationships  between  different  parts 
of  the  stream  in  Bozeman  Creek  (Figures  4  and  5). 

In  Bozeman  Creek,  number/km  for  station  36,  a  control  station  of  greater- 
than-average  length  and  lying  just  below  a  mink  farm  where  mink  feed  and  manure 
are  reportedly  introduced  into  the  stream,  was  significantly  higher  than  in 
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Table  5.   Estimated  number  of  over-20-cm  trout/km  for  individual  stations. 


Urban 

Number  trout 

or 

over  200  mm 

95%  Confidence 

Stream 

Stations 

Control 

per  km 

Interval 

Bozeman  Creek 

7 

C 

644 

469  - 

1030 

8 

C 

457 

283  - 

914 

9 

C 

533 

371  - 

857 

13 

C 

611 

506  - 

802 

14 

C 

412 

266  - 

714 

15 

C 

542 

374  - 

869 

19 

U 

209 

146  - 

428 

20 

U 

499 

297  - 

940 

21 

u 

386 

193  - 

914 

22 

U 

272 

212  - 

363 

23 

U 

549 

435  - 

795 

24 

U 

240 

200  - 

400 

26 

U 

319 

150  - 

798 

27 

U 

296 

245  - 

459 

28 

U 

135 

127  - 

195 

36 

C 

882 

729  - 

1128 

60 

C 

326 

289  - 

400 

61 

C 

299 

269  - 

359 

Mathew  Bird  Creek 

1 

u 

253 

» 

2 

u 

167 

130  - 

315 

3 

u 

179 

« 

4 

u 

155 

141  - 

247 

5 

c 

30 

» 

6 

C 

121 

111  - 

192 

7 

c 

87 

» 

Figgins  Creek 

1 

u 

64 

43  - 

133 

2 

C 

241 

204  - 

320 

Fleshman  Creek 

1 

u 

70 

» 

2 

u 

105 

98  - 

173 

3 

C 

83 

63  - 

136 

Not  available 
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Table  6.   Numbers  of  over-20-cm  trout/km  for  grouped  stations  (sections). 


Urban 

Number  trout 

or 

over  200  mm 

95%  Confidence 

Stream 

Stations 

Control 

per  km 

Interval 

Bozeman  Creek 

7 

-  9 

C 

579 

452  - 

757 

13 

-  15 

C 

511 

433  - 

726 

19 

-  21 

u 

355 

251  - 

570 

22 

-  24 

u 

351 

305  - 

442 

26 

-  28 

u 

218 

190  - 

294 

36 

c 

883 

729  - 

1128 

60 

-  61 

c 

311 

284  - 

352 

Mathew  Bird  Creek 

1 

-  3 

u 

183 

143  - 

287 

4 

u 

155 

141  - 

247 

5 

-  7 

c 

79 

76  - 

96 

Figgins  Creek 

1 

u 

64 

43  - 

133 

2 

c 

241 

204  - 

320 

Fleshman  Creek 

1 

u 

70 

Numbers 

to  low 

to  compute  CI 

2 

u 

105 

98  - 

173 

3 

c 

83 

63  - 

136 
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any  of  the  grouped-station  sections  except  section  7-9,  which  itself  had 
significantly  greater  trout  abundance  than  several  of  the  other  sections  (Figure  5). 

The  relative  abundance  of  rainbow  trout  by  number  decreased  in  the  upstream 
direction,  from  about  80%  of  the  trout  population  in  the  lowermost  section  to 
about  49%  in  the  uppermost  section  (Table  7).   Brook  trout  rose  in  relative 
abundance  further  upstream,  their  highest  proportion  being  51.4%  of  the  trout 
population  in  section  60-61. 

Few  brown  trout  were  captured  in  Bozeman  Creek.   Station  60  yielded  the 
most — four  with  total  weight  of  0.673  kg.   Station  9  had  the  largest-sized 
brown  trout — three  for  a  total  weight  of  1.374  kg.   Mountain,  longnose,  and 
white  suckers  were  collected  as  a  large  proportion  of  the  fish  populations 
in  stations  7  through  14  and  21  through  24. 

Brown  trout  was  the  principal  species  collected  in  the  urban  sections 
of  Fleshman  Creek.   Other  species  captured  in  urban  stations  were  rainbow 
and  brook  trout.   Brook  trout  (population  data  from  Clancey  and  Gould 
1982)  was  the  only  species  collected  in  the  Fleshman  Creek  control  station. 
Comparisons  with  Other  Streams 

Trout  abundance  varies  greatly  throughout  small  streams  of  the  Gallatin, 
Madison,  and  Jefferson  River  drainages.   There  are  34  population  estimates 
available  on  31  rather  scattered  stream  sections  with  lengths  of  about  130 
to  380  m  on  17  creeks  (Table  8).   Trout  biomass  per  unit  length  of  stream 
averaged  66  kg/km  and  ranged  from  8.9  to  222  kg/km,  with  50%  of  cases  lying 
between  27  and  87  kg/km.   With  those  sections  omitted  which  are  known  to  have 
been  physically  altered  by  artificial  means  (straightening  or  bank  landfill, 
urban  and  otherwise),  the  mean  of  the  20  more  natural  remaining  cases  was  77.5 
kg/km,  and  50%  of  the  cases  lay  between  42  and  112  kg/km. 
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Table  7.   Species  composition  of  the  trout  pouplations  in  streams  of 
the  Bozeman  Creek  system,  August  1982. 


Rainbow  Brook 

Section  trout  trout  Total 


Bozeman  Creek 

7-9               200  54                 254 

78.7%  21.3% 

13  -  15              170  46               216 

78.7%  21.3% 

19  -  21               121  20                 im 

86  %  14  % 

22  -  24              137  35                172 

79.7%  20.3% 

26  -  28               62  39                101 

61.4%  38.6% 

36                88  48               136 

64.7%  35.3% 

60                35  47                82 

32  %  68  % 

61                 68  62               130 

52.3%  47.7% 

60  -  61              103  109                212 


1-3               30  47                77 

39  %  61  % 

4                0  43                43 

0  %  100  % 

5-7                4  25                29 


1  0  26  26 

2  0  70  70 


200 

54 

78.7% 

21.3% 

170 

46 

78.7% 

21 . 3% 

121 

20 

86  % 

14  % 

137 

35 

79.7% 

20.3% 

62 

39 

61.4% 

38.6% 

88 

48 

64.7% 

35.3% 

35 

47 

32  % 

68  % 

68 

62 

52.3% 

47.7% 

103 

109 

48.6% 

51.4% 

Mathew 

Bird  Creek 

30 

47 

39  % 

61  % 

0 

43 

0  % 

100  % 

4 

25 

13.8% 

86.2% 

Figgins  Creek 
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Table 


Trout  abundance  and  stream  characteristics  in  various  small  streams  of 

the  Gallatin,  Madison,  and  Jefferson  River  drainages,  southwestern  Montana. 


Mean 

Drain- 

Yield 

Mean 

Sta- 

Distance 

Trout 

Alt- 

discharge 

age 

(liters 

bed 

ion 

above 

standing 

era- 

at  mouth 

area 

per  sec 

slope 

1th 

mouth 

crop 

Data 

Stream      tion* 

(mVsec) 

(km^) 

per  km^) 

(m/km) 

(m) 

(km) 

(kg /km) 

Sourc 

Bozeman  Cr. 

1.02 

168 

6.1 

19 

MDFWP  sta 

24 

305 

13.7 

41.7 

1 

Sta  60-61 

10 

363 

6.5 

60.4 

2 

Sta  36 

10 

130 

3.9 

155.0 

2 

Sta  26-28  U 

10 

316 

3.0 

42.6 

2 

Sta  22-2A  U 

10 

305 

2.3 

75.0 

2 

Sta  19-21   U 

10 

298 

2.0 

53.8 

2 

Sta  13-15   S 

10 

321 

1.4 

113.2 

2 

Sta  7-9 

10 

330 

0.8 

108.9 

2 

Mathew  Bird  Cr. 

0.09** 

? 

? 

33 

Sta  1-3    U 

300 

0.6 

35.4 

2 

Sta  4      U 

142 

1.5 

35.7 

2 

Sta  5-7    ? 

302 

3.0 

24.7 

2 

Figgins  Cr. 

0.001** 

7 

7 

13.5 

Sta  1      U 

188 

0.2 

15.3 

2 

Sta  2 

191 

0.5 

51.0 

2 

Fleshman  Cr. 

0.07** 

7 

9 

23 

Sta  1      U 

157 

ca  1 

23.0 

2 

Sta  2      U 

133 

ca  2 

17.3 

2 

Sta  3      S? 

300 

ca  3 

18.6 

3 

Hell  Roaring  Cr. 

1.27 

78 

16.3 

65 

305 

ca  1 

26.8 

1 

Hyalite  Cr. 

1.97 

306 

6.4 

29 

305 

ca50 

111.6 

1 

Porcupine  Cr. 

0.69-0.86 

78 

9-11 

51 

305 

ca  1 

8.9 

1 

S,  Cottonwood  Cr. 

1.01 

109 

9.3 

41 

305 

ca  14 

163.7 

1 

S.  Fk.  Spanish  Cr 

7 

105 

? 

60 

305 

ca  12 

26.8 

1 

Squaw  Cr. 

1.37 

142 

9.6 

53 

305 

ca  4 

66.9 

1 

Taylor  Fork 

2.77-3.01 

277 

10-11 

27 

305 

ca  0.5 

61.0 

1 

W.  Fk.  Gallatin  F 

.   1.95 

202 

9.6 

23 

305 

ca  1 

26.8 

1 

Rocky  Cr. 

1.01 

225 

4.5 

9.5 

Sect  A 

282 

4.8 

193-222 

4 

Sect  B     S 

381 

1.6 

25-87 

4 

Sect  C     S 

305 

0.8 

72-82 

4 

South  Boulder  R. 

1.45 

246 

5.9 

28 

305 

7 

59.5 

1 

Whitetail  Ci . 

0.62 

482 

1.3 

22 

305 

7 

159.1 

1 

Jack  Cr. 

1.00 

166 

6.0 

36 

305 

? 

47.6 

1 

N.  Meadow  Cr. 

1.02 

137 

7.4 

39 

305 

7 

35.7 

1 

*U  =  urban,  S  =  straightened 

Data 

sources: 

1 = MDFWP  1981 

3  =  Clancy  &  Gould 

1982 

**Summer  low  disch 

arge 

2  =This 

study 

4=  Wells 

1977 
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Within  a  single  stream,  Rocky  Creek  near  Bozeman,  there  was  variation 
from  25  kg/km  in  a  straightened  section  in  springtime  to  222  kg/km  in  an 
unaltered  section  in  autumn.   In  Bozeman  Creek,  where  we  have  data  from  8 
sections  (encompassing  18  shortered  stations),  trout  abundance  ranged  from 
42  kg/ km  in  one  urban  section  and  in  one  natural  section  far  upstream  to 
155  kg/km  in  a  relatively  unaltered  section  (station  36)  which,  however, 
may  have  been  influenced  by  nutrients  from  a  mink  farm. 

Trout  biomass  abundance  of  natural  stream  sections  was  inversely  related 
to  channel  steepness  (Table  8,  Figure  6).  From  the  plot  of  trout  abundance 
against  channel  gradient  in  Figure  6,  it  is  evident  that  urban  or  straightened 
sections  held  less  trout  than  natural  stream  sections  of  equivalent  gradient. 
The  estimates  of  stream  gradient  must  be  regarded  as  rather  rough.   Within 
the  apparent  gradient-trout  relationship,  sections  of  Bozeman  Creek  having 
highest  abundance  of  trout  compared  favorably  with  abundances  in  other 
streams  that  had  natural  channels . 
Further  Analysis  Planned 

Further  refinement  and  augmentation  of  various  aspects  of  the  data 
analysis  in  this  study  are  warranted.   Confidence  limits  on  estimates  of 
biomass  remain  to  be  calculated,  as  well  as  improvements  in  computations 
involving  grouped  data  for  contiguous  stations  (stream  sections).   HQI 
analysis  and  multiple-correlation  analysis  of  data  from  the  30  stations 
should  be  repeated  using  stream  section  data  instead.   Comparisons  of 
means  should  be  further  analysed  with  weighting  of  data  according  to 
stream  station  or  section  length. 

Most  aspects  of  trout  population  size  structure  and  age-growth  analyses 
remain  to  be  performed.   Scale  samples  were  taken  from  most  fish  captured 
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Figure  6.   Relationship  of  standing  crops  of  trout  to  stream 
gradient  in  creeks  of  the  Gallatin,  Madison  and 
Jefferson  River  drainages  in  southwestern  Montana 
(data  from  Table  8) . 
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during  electrofishing,  but  these  have  not  yet  been  mounted  and  read. 

It  is  expected  that  the  further  analyses  will  result  primarily  in 
augmentation  of  present  information,  not  in  revision  of  the  primary 
conclusions. 
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Appendix  I.   St 

.udy  station  marker  locat 

ions  (upstream  bou 

ndaries  of 

Stat 

ions) . 

Landmark , 

Street  or 

Property 

Station 

Designation 

Location 

Bozeman 

Creek 

7 

Bond  St. 

NEi, 

SWi, 

NEi, 

T2S, 

R6E, 

S6 

8 

Midway 

SEi, 

swi, 

NEi, 

T2S, 

R6E, 

S6 

9 

Gold  St. 

NEJ, 

NWi, 

SEJ, 

T2S, 

R6E, 

S6 

13 

NWi, 

NWi, 

SEJ, 

T2S, 

R6E, 

S6 

14 

Kenyon-Noble 

swi, 

NWJ, 

SEi, 

T2S, 

R6E, 

S6 

15 

Tamarack 

swi, 

NWi, 

SEi, 

T2S, 

R6E, 

S6 

19 

Kwik-Way 

swi, 

swi, 

SEi, 

T2S, 

R6E, 

S6 

20 

NWi, 

NWi, 

NEi, 

T2S, 

R6E, 

S7 

21 

Rouse 

NWJ, 

NW^, 

NEi, 

T2S, 

R6E, 

S7 

22 

Skiway 

SEi, 

NEi, 

NWi, 

T2S, 

R6E, 

S7 

23 

NEi, 

SEi, 

NWi, 

T2S, 

R6E, 

S7 

24 

Library 

NEJ, 

SEi, 

NW4, 

T2S, 

R6E, 

S7 

26 

Pavillion 

NEi, 

NE4, 

swi, 

T2S, 

R6E, 

S7 

27 

Bogert 

SE4, 

NE4, 

swi, 

T2S, 

R6E, 

S7 

28 

Story 

SEi, 

NEi, 

SWi, 

T2S, 

R6E, 

S7 

36 

Ice  Pond 

swj, 

NWi, 

NEi, 

T2S, 

R6E, 

SIB 

60 

Picton 

NWJ, 

SEh, 

NEi, 

T2S, 

R6E, 

SIB 

61 

Meadowlark 

swi, 

SE4, 

NEi, 

T2S, 

R6E, 

S19 

Mathew  Bird  Creek 


Teepee 


Wood  Brook 
Graf  1 
Graf  2 
Graf  3 


NWJ,  NEi,  NWi,  T2S,  R6E,  SlB 

NWi,  NEi,  NWi,  T2S,  R6E,  SlB 

swi,  NEi,  NWi,  T2S,  R6E,  SlB 

NE^,  NWi,  swi,  T2S,  R6E,  SlB 

NEi,  SWi,  NWi,  T2S,  R6E,  S19 
El,    swi,  NWi,  T2S,  R6E,  S19 

SEi,  swi,  NWi,  T2S,  R6E,  S19 


Figgins  Creek 


Hoffman 
Control 


SEi,  NWi,  swi,  T2S,  R6E,  SlB 
NWi,  SWi,  SWi,  T2S,  R6E,  SlB 


Fleshman  Creek 


Sacajawea 

Crawford 

Clancy 


NWi,  SWi,  NWi,  T2S,  R9E,  S24 
SEi,  NWi,  NWi,  T2S,  R9E,  S24 
SEi,  SEi,  swi,  T2S,  R9E,  Sl4 
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Appendix  III.  Multiple-regression  models  describing  numerical  density  of  over-20-cm  trout  (trout/km) 
as  a  function  of  stream  nabitat  variables. 

1 .  All  30  stations  included  (urban  and  controls) 

log  trout/km  =  15.289  +  1.77  log  water  velocity 

+  0.153  log  critical  period  flow/mean  flow 

+  0.317  log  of  solid  overtiead  cover 

-  0.3'»8  log  max  summer  water  temp 

+  0.2't3  log  vegetational  substrate  rating 

r  =  0.8753,  simple  r^  -  0.763,  adjusted  r^*  =  0.713 

Order 

or  Adj.  r     F-values  at  step 

entry     Variable      at  step    log  wtrvel     log  CPf/MF log  cover log  maxT    log  substi 


log  wtrvel 
log  CPF/MF 
log  cover 
log  maxT 
log  substr. 


.5't2 

.670 
.703 
.713 


35.297 
23.606 
25.99't 
30,856 
25.552 


8.768 
8.6'.7 
If. 902 
1.15't 


5.352 

6.12't 
5.105 


3.8't5 
5.902 


1.922 


2.  Urban  Stations  Only 


log  trout/km  =  -  64.255  -  0.762  log  water  velocity 

+  0.905  log  vegetational  substrate  rating 

+  1.560  log  NO  -N 

+  0.72't  log  annual  streamflow  variation 

+  1.827  log  max  T 

+  l.O'i't  log  critical  period  flow/mean  flow 

-  0.'«22  log  %  pool  &  turbulence  cover 

+  0.215  log  %   eroding  bank 


r  -   0.976,  simple  r 

Order 
or 
entry  Variabl 


0.952,  adjusted  r 


Adj.  r 


F-value  at  step 


at  step  log  wtrvel  log  substr  log  NO   log  ASFV  log  maxT  log  CPF/MF  log  pool  loger.  bu- 


1    log  wtrvel 

59't 

22.966 

, 

2    log  substr 

710 

35.277 

6.567 

3    log  NO 
if    log  ASrv 

715 

28.5'.0 

7.687 

1.252 

7't6 

IZ.'tO't 

5.598 

5.7'.7 

2.'.66 

1 

5    log  maxT 

7'tl 

3.051 

5.828 

'..509 

2. 1.60 

0.796 

6    log  CPF/MF 

8'*5 

0.872 

15.550 

5.'t63 

3.519 

8.785 

7.'.7'. 

7    log  pool 

858 

0.829 

19.159 

7.55'. 

1.961 

10.756 

't.567 

1.967 

8    log  er.buk. 

898 

3.5't3 

50.580 

1'..650 

2.778 

19.088 

7.7'.0 

5.225 

-..litS 

*  Adjusted  R  r  r^  - 

[(K- 

1)/N 

-K)][l-r2] 

,  where  K  := 

number  of 

independent 

variables 

N  =  number 

of  cases. 
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Appendix  III.  Continued. 


3.  Control  Stations  Only 

log  trout/km  =  hk.Zll   +  1.300  log  water  velocity 

-  0.768  log  maxT 

-  0.3't't  log  annual  streamflow  variation 
+  0.298  log  NO  -N 

r  =  O.g'tg,  simple  r^  =  0.901,  adjusted  r^  =  0.857 


Adj.  r^ 
at  step 

.578 
.761 
.837 
.857 

F-value  at  step 


log  wtrvel       log  maxT 


log  ASFV      log  NO 


1  log  wtrvel 

2  log  maxT 

3  log  ASFV 
h  log  NO 


40.700 
61.105 
45.61'. 


10.162 
20.625 
20.266 


6.114 
7.032 


2.385 


I 


